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We analyze charge tunnehng statistics and current noise in a superconducting single-electron 
transistor in a regime where the Josephson-quasiparticle cycle is the dominant mechanism of trans- 
port. Due to the interplay between Coulomb blockade and Josephson coherence, the probability 
distribution for tunneling events strongly deviates from a Poissonian and displays a pronounced 
even-odd asymmetry in the number of transmitted charges. The interplay between charging and 
coherence is reflected also in the zero-frequency current noise which is significantly quenched when 
the quasi-particle tunneling rates are comparable to the coherent Cooper-pair oscillation frequency. 
Furthermore the finite frequency spectrum shows a strong enhancement near the resonant transition 
frequency for Josephson tunneling. 

PACS numbers: 72.70.+m, 73.23Hk, 74.50+r 



I. INTRODUCTION 

Shot noise in a mesoscopic conductor is a consequence 
of the stochastic character of electron tunneling and of 
the discreteness of charge. Unlike thermal noise, shot 
noise describes the non-equilibrium fluctuations of cur- 
rent; therefore, the study of current fluctuations can pro- 
vide further understanding of properties related to cor- 
relation mechanisms, internal energy scales or the carrier 
statistics whip|h^eannot be obtained by measuring the av- 
erage currenltl'ElH'a. 

A well studied example of physical processes where 
electron correlations play a dominant role is the phe- 
nomenon of Coulomb blockade. In a system of small 
tunneling junctions, due to the large electrostatic energy 
(as compared to temperature or voltages), the electronic 
charge is transported one by one. This effect leads to 
many remarkable features in transport properties and ] 
been a subject of extensive study for the last decades 
As an example, the strong dependence of the current- 
voltage characteristics on the gate charge was exploited 
to use a single electron ^transistor (SET) as a highly 
sensitive charge detectorQ, and proposed as measuring 
apuar|atpj of the charge state of a Josephson quantum 
bit@ot2ltll. Since it leads to a strong correlation of con- 
secutive tunneling events. Coulomb blockade has turned 
out to manifest itself in a peculiar way on the current 
noise. Such aixpjfiEQcLhas been studied both in the secj 
tial tunnelingBDLSO and in the cotunneling regime 
Additional interest in studying noise in single electron 
devices comes from the recently proposed schemes that 
employ current fluctuation measpjements to detect en- 
tanglement in solid state systemst^. 

An even richer scenario occurs when the coherence of 
charge carriers is maintained over a significant portions 



of the system. Such a circumstance is encountered quite 
often when the tunneling junctions are superconducting. 
In this case, the charge carriers are Cooper pairs, and 
their coherent tunneling across the junctions gives rise to 
a series of pronounced stjij4ij:itj^ in the I-V characteris- 
tics at sub-gap voltageglZHl3't3. Furthermore, as analyzed 
in Ref. |2^, the scattering of quasi-particles (and conse- 
quently the shot noise) in superconducting point contacts 
is significantly enhanced in the presence of the supercur- 
rent produced by a coherent flow of Cooper pairs. 

In this paper, we analyze a superconducting double 
tunnel junction device, operating in a suitably chosen 
bias voltage regime, such that one of the junctions of 
the SET is on resonance for Cooper pair tunneling (the 
case where Cooper pair resonance occurs on both junc- 
tions has been recently analyzed in Ref. ^l|). The inter- 
play between coherence and interaction is explored by 
sweeping the operating point of the device through the 
Cooper pair resonance. We will show that the fluctua- 
tions of the charge on the central island are sensitive to 
both Coulomb blockade and quantum coherence. More 
pronounced effects arise in the regime in which the rates 
of incoherent quasi-particle tunneling matches the fre- 
quency of coherent Cooper-pair oscillation. This gives 
rise to an enhanced fluctuation of charge in the central- 
island and to a substantial suppression of the current 
noise. By investigating the statistics of the tunneling 
events, we show that the suppression in the shot noise is 
related to the deviation of the counting statistics from 
the Poissonian distribution. The probability distribu- 
tion of tunneling events exhibits a parity dependence and 
remains non-Poissonian in a wide range of parameter 
values. The interplay between coherence and Coulomb 
blockade affects the overall charge transport and is also 
clearly observed in the finite frequency behavior of the 
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FIG. 1: Schematic diagrams of (a) the superconducting 
SET device and (b) the transition processes between relevant 
charge states. 



current noise. Its power spectrum displays a sharp res- 
onance peak at the Josephson frequency, resulting from 
coherent oscillations between two quantum states. 

The work presented here applies to the setup used in 
a recent experimentE3 to probe the coherent evolution 
of quantum states irt-a Cooper pair box as well as in 
an earlier experimentllll on resonant Cooper pair tunnel- 
ing. In this paper we extend the results of Ref. |2^. The 
paper is organized as follows. In Section II, we intro- 
duce the model to describe the SET transistor and we 
describe the relevant processes involved in the Josephson- 
quasiparticle cycle. In the same section we also introduce 
the master equation, in its general form, used to obtain 
all the results of this work. The solution of the master 
equation is worked out to analyze in detail the fluctua- 
tions of the charge on the central electrode and, in partic- 
ular, the spectrum of the charge fluctuation as measured 
by a detector coupled to the island (Section |l|), the 



counting statistics for the tunneled charges (Section |IV|), 
and the shot noise of the current (Section ^). In Sec- 
tion 1^ we discuss in some details the properties of the 
shot noise at zero frequency, making explicit the results 
already contained in the counting statistics. In the same 
section we also discuss the frequency dependence of the 
shot noise. Section VI is devoted to the conclusions. 



II. THE MODEL 

The system we consider is a superconducting SET (see 
Fig. [^), which consists in a small central electrode (is- 
land) connected by tunnel junctions to two leads and 
capacitively coupled to a gate electrode. The electro- 
static energy of the island can be adjusted by control- 
ling the gate voltage Vg. A transport voltage V is ap- 
plied to the outer leads, determining the current flowing 
through the device. Letting Cq^Cl and Cr be the gate 
and left and right junction capacitances, respectively, the 
electrostatic charging energy is given by Ec = I^Cy. 
where Cy. = + Cl + Cr is the total capacitance of 
the island. The device operates in the Coulomb blockade 
regime, i.e. with the charging energy much larger than 



both the Josephson coupling energy Ej and the thermal 
energy ksT {Ec ^ Ej^ksT)^ but still much smaller 
than the superconducting gap A which we suppose to be 
the largest energy scale in the problem. Charge can pass 
through the tunnel barriers coherently and incoherently. 
In addition to quasi-particle tunneling, Josephson effect 
allows to maintain coherence in the Cooper-pair tunnel- 
ing processes. This coherence does play an important 
role in transport as long as charge states which differ by 
one Cooper pair are almost degenerate. The Hamiltonian 
of the system is given by 



Hl^Hr^Hi^Ht^ He 



(1) 



where i^Q, {a = L, I) is the BCS Hamiltonian of the 
left (L), right (R) leadj-guid of the central island (I). The 
tunneling HamiltonianEj is 



j=L,R kqa 



kq^ 



h.c. 



(2) 



where T^q is the tunneling amplitude and c^^ (cka) cre- 
ates (annihilates) a particle with momentum k and en- 
ergy Ska in electrode a. The variable ^L,i? is the super- 
conducting phase difference at the left (right) junction 
and it is canonically conjugated to the number til^r of 
electrons that have passed across the left (right) junc- 
tion out of the central electrode ([n/c,^j] = iSjk). Fi- 
nally, He is the electrostatic energy. 



He = Ec{n^no)^ ^eVuR. 



(3) 



where n = —ul — tlr is the number of excess electrons on 
the central island, while eno = CrV -\- CgVg is the offset 
charge due to the applied voltages. 

For later convenience it is useful to define the part of 
the Hamiltonian which accounts for the coherent dynam- 
ics of the macroscopic variable n. It includes the charging 
and the Josephson terms. By properly adjusting the bias 
and gate voltages, one can put either the right or left 
junction at such a resonance for Cooper-pair tunneling. 
We consider the case of resonance across the left junc- 
tion, and consequently we keep only the corresponding 
Josephson coupling 

Ho = Ec{n + nof + eVuR - Ej cos(2^l) • (4) 

A quasi-particle tunneling into {out of) the island 
across the junctions leads to the transition n ^ n + 1 
(n n — 1) of the island charge. The rates of these 
incoherent processes are given by the relation 



- L/R 



(n) = coth(p^^^^^ ) ± 1 



■n,± ) 



2e 



(5) 



where S^^^ 

Ec {m — n){m -\- n -\- 2no) , and Iqp{E) is the quasi-particle 
tunneling current at the bias volta ge E/e in the absence 
of charging effects (see, e.g., Ref. 
pression of /qp in terms of Tj^q). 



25 for an explicit ex- 
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Since the SET transistor operates in the charge regime 
{Ec ^ Ej)^ we can take advantage of the strong sup- 
pression of charge fluctuations to use the eigenstates of 
n as basis states for the island. Moreover, we focus on 
the bias regime \eV\ 2 A + Ec where only the two 
charge states with n = and n = 2, are nearly degener- 
ate. Such a condition implies that quasi-particle tunnel- 
ing only takes place from the central island toward the 
right electrode, while the left junction allows only for co- 
herent Cooper pair tunneling. Furthermore, we suppose 
that the Josephson energy of the right junction is negli- 
gible (the corresponding term has already been omitted 
from iJo, which is justified within the rotating wave ap- 
proximation). All these conditipas are met in the recent 
experiment by Nakamura et a/.EHl, designed to probe the 
state of the island via the detection of the incoherent 
tunneling current. In this situation one can imagine that 
the coherent Cooper-pair tunneling occurring across the 
left junction is interrupted from time to time by quasi- 
particle tunneling across the right junction, as sketched 
in Fig.0 (b). 

Due to the strong Coulomb blockade, it suffices to keep 
the three charge states, n = 0,1,2, and two tunneling 
rates, Fi = F^(l) and F2 = F^(2); the other tunneling 



rates are exponentially suppressed. In order to simplify 
the notation, we assume that Fi = F2 = F, which is a 
very good approximation in the regime we are interested 
in. For example, in the experiment of Ref. |2^, 1/Fi = 8ns 
and I/F2 = 6ns. 

The transport properties of the system in the set-up 
described above can be well described in terms of two 
variables, either n and ul or n and ur {n = —ul — tlr). 
However, the quantum dynamics of these system is af- 
fected by the quantum noise due to the coupling to the 
environment provided by the fermionic bath. In order 
to describe this effect, we adopt a master equation ap- 
proach, which4|ias been widely used to describe quantum 
open systemsE^. A master equation for the reduced den- 
sity matrix p{t) = TrqpPtot(^) is obtained by taking the 
trace over the fermionic degrees of freedom from the Li- 
ouville equation {% = 1) 

dtptot (t) = -i [Htot , ptot (t)] (6) 

for the density matrix ptot of the system plus environ- 
ment. The resultmg emiation can then be written in the 
Lindblad form astSLBE^'EjO 



dtp{t) 



[Ho, Pit)] + :^ E [2i„ Pit) Li - LlL„ p{t) - p{t) it L„ 



= 1,2 



(7) 



Here Ln is a Lindblad operator corresponding to the 
quantum jump n ^ n — 1 and ur n^^ + 1, i.e., in 
the |n, ni^)-basis = \n — 1, ur + 1) (n, nR\. The first 
term describes a purely phase-coherent dynamics, while 
the second one is responsible for both dephasing and re- 
laxation due to the quasi-particle tunneling. 

The solution of the Eq. behaves in distinct ways in 
the two limiting cases of strong and weak coupling with 
the quasi-particle reservoir. In the strong dephasing limit 
(either V^Ejoie^Ej^ see below), the dephasing 
time r^^, which describes the decay of the off-diagonal 
elements of p to their stationary values, is small compared 
to the relaxation time which sets the time-scale for the 
variation of the diagonal elements (i.e. population of the 
charge states). The relaxation time is given by 



2E]T 



(8) 



On the other hand, in the weak dephasing limit (F, £ <C 
£'j), there is no such a clear separation of time scales; 
both the diagonal and off-diagonal elements vary over 
the same time scale l/F. 



III. FLUCTUATIONS OF THE CHARGE ON 
THE ISLAND 

A first insight into the interplay among coherent 
Cooper-pair tunneling. Coulomb blockade, and incoher- 
ent quasi-particle tunneling can be obtained by examin- 
ing the fluctuations of the charges on the island as a func- 
tion of the quasi-particle tunneling rate F, the gate volt- 
age, and the Josephson coupling energy. In this case, we 
only need to keep track of the variable n, and thus define 
a reduced density matrix for the central island charge, 
crit) = Tin^ p{t)^ which satisfies an equation identical to 
Eq. (|^), but with Ln =T-\p — 1) (^| now operating on 
the reduced 'n' space onlyEZl. In the stationary state, the 
master equation has the solution {(Jmn = {m\a\n)) 



i + (4£2 + r2)/E3 



3 + (4^2 



CTU (722 



(^02 = 0^20 



1 



3 + (4£2 + T^)/E] 
_ . Ej{(Joo — (J22) 
~* T + 2is 



(9a) 
(9b) 
(9c) 



Here s = 4£^c'(l + no) measures the energy difference 
between the state with n = and n = 2 charge on the 



.J I IJ 2 2.5 



.5 S i.5 2 i-J 



FIG. 2: Variance ((An)^) of the charge on the island as a 



function of Fi 



F and e. 



island; the Cooper pair resonance corresponds to e = 0. 

From the stationary-state solution (||), we evaluate the 
characteristic function C(^) = (^e~^^^ for the quantum 
variable n. It is given by 



c{0) = 



1) + + r2 



4^2 + r2 + sEj 



(10) 



From C{0) one can evaluate all the statistical moments 
of the charge on the island; in particular, we concentrate 
on the variance 



((An)2) 



5(4e2 



■6Ej 



[4£2 + r2 



(11) 



shown in Fig. 

An interesting point about the result in Eq. (|ll] ) is that 
the maximum of the variance is found at Fopt — \/3/5 Ej 
for e = 0. In other words, the fluctuation is enhanced 
when the decoherence time matches the time for the 
Josephson coherent oscillations. Moreover as |5| increases 
(up to l^l < Ej y^3/20), the optimal value of F decreases 
as Fopt — \/3£^j/5 — 46^; meaning that e enhances the 
effective dephasing rates. These two features will appear 
more clearly when we discuss the statistics (Section IV) 
and noise (Section ^ of the transport across the junc- 
tions. 

An important quantity to consider is the fluctuation 
spectrum for the number of electron charge residing on 
the island. As already discussed by many authors, it 
is important kj^^jmclude the back-action of the measur- 
ing apparatud§l3ll3, which could be a SET transistor ca- 
pacitively coupled to the central island. We include the 
charge detector coupling via an Hamiltonian term of the 
form 



det 



-fiD 



where D is a detector operator. Assuming the correlation 
time for the detector to be the fastest time scale of the. 
problem, we write (here we follow Averin's treatmentE£l 



{D{t + T)D{t))=-id5{T) 



(12) 



The non-zero value of is the essential cause of the 
measurement back-action. Indeed, a term proportional 



FIG. 3: Charge fluctuation spectrum in the weak dephasing 
regime at degeneracy (s = 0), with the indicated back-action 



rate values and with T/Ej 
(right). 



10-^ (left) and V/Ej = 0.1 



to enters the master equation Eq.(|^), thus affecting 
the time evolution of the system variables. To effectively 
measure the charge number, we look at an output detec- 
tor operator O, which, in the linear regime, evolves as 
0{t) = 0^°) + An(t), (see Ref.|l^). The response coeffi- 
cient A is determined by the imaginary part of the equilib- 
rium correlation function {0{t -\- T)D(t)) . Furthermore, 
A can be related to jd so that we can write for the signal 
to noise ratio 



1^2jdSn{u;) 



(13) 



where it is assumed that the real part of the O-D cor- 
relator vanishes (which is the most favorable case for a 
measurement). Here Sn{^) is the charge number fluctu- 
ation spectrum evaluated as 



e'^^^dr (14) 



where the time evolution is obtained from a modified 
master equation including the back-action. Note that 
here we use the symmetric correlation function since the 
island charge itself is coupled to the quasi-particle bath, 
so that the detector can also receive energy from the sys- 
tem. 

As shown in Fig. (^, the spectrum displays a reso- 
nance peak at the Josephson frequency which is broad- 
ened by the quasi-particle rate F (the peak is only visible 
in the weak dephasing regime, otherwise it is completely 
washed out independently of the value of 7^^). As the 
back-action rate jd increases, a maximum develops at 
zero frequency, which finally hides the resonance struc- 
ture. This enhanced zero frequency noise results from 
incoherent transition induced by the detector coupling. 

Near the maximum at :^ Ej, and for F, 7^^ <C Ej^ 
the spectrum takes the approximate form 



-Ef 



(r + 27d) 



(15) 
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IV. COUNTING STATISTICS FOR THE 
TRANSMITTED CHARGE 

We turn now to a description of the statistical distri- 
bution of the number of clwges transmitted through the 
system during a period r,E3. Specificahy, we wih exam- 
ine the probabihty Pt{N^T) that N electrons have been 
transferred across the right junction during the interval 
[t, t + r]. We note that 

Pt{N, r) = Y,p{N + ni^, t + r; ur, t) , (16) 

riR 



where p{nR -\- N^t-\-r] ur, t) is the joint probability that 
riR electrons have passed across the right junction up to 
the time t and ur -\- N electrons up to time t -\- r. To 
obtain Pt (A/", r) , we define a characteristic matrix 



Gt{0,r) = y exp(-z^iV)Tr(ni^ + 7V|e-^^-^n^i?)(^i?l/>tote+^^^°^n^i? + ^) 



(17) 



defined so that Tr Gt(6>,r) 



AdN\ 



is the characteristic function for Pt{N^r). Namely, 



PtiN,T) 



r — 



,+i9N 



TvGt{0,T), 



where the trace is taken over the states |n). Following the same procedure that led to Eq. 
Gt{0,T) satisfies the following master equation: 



drGt = -i [Ho, Gt] + l Yl [^e'^L^GtLi - GtllL^ - liL^Gt] 



(18) 

one can show that 
(19) 



= 1,2 



with the initial condition Gt{0,0) = p{t) \nR). 

Here we will consider two limiting cases for the solution, the strong and the weak dephasing limit (see the discussion 
at the end of Section II). We find that in the strong dephasing case (T ^ Ej or s ^ Ej) 



TTGt{e,T) [aoo{t) + zan{t) + z^a22{t)] exp 
while in the weak dephasing limit at resonance (T Ej, e = 0) 



(1-.2) 



(20) 



TTGt{e,T) = 



cosh ^-iM + l + 2an(^)+2.[l-.n(^)] ^.^^ Fr^ 



+ — {I- z)lmao2{t) 
Ej 



f{z) 



cosh^:IM+li±i3inh^:IM 



f{z) 



■ exp 



2,Tt \ 



{1- z)< Im(J02(t) COs(£;jt) - - [(Joo(t) - (T22{t)] sm{EjT) 



I exp ^- 



(21) 



where z = e*^, and f{z) = Vl + 8z^. In Figs. || and ||, dephasing limit for various starting conditions cr(t). In 
the resulting statistics are shown for the weak and strong the limit of strong dephasing, the counting statistics de- 
dephasing cases, respectively, in the transient state (i.e. pends sensitively on the initial state cr(t). This has been 



t 0). 



exploited in Ref. 22, where the measurement of the quan- 



As seen in Eqs. ^) and (H), Gt{0,r) and hence tum state is performed with the system taken far from de- 

Pt{N, r) depend on the charge state of the island at time generacy. On the contrary, the dependence on the initial 

t. This point is further illustrated in Fig. | where the condition is quickly lost in the weak dephasing regime, 

counting statistics at short time is shown in the strong since the strong Josephson energy can rapidly produce a 
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P(N) 




FIG. 4: Pt^o(N,r) obtained by numerically solving equation 
(^ij), in the case of Fr = O.IF, with initial condition poo(O) = 

1,P22(0) =0. 




FIG. 5: The short time probability for the first few tunneling 
events in the weak dephasing regime (F = O.lEj) and at res- 
onance (e = 0). In this regime, the counting probability does 
not depend on the initial condition for p 



change in the state, before quasi-particles have any time 
to be produced. 

Another important limit to consider is the stationary 
state {t oo), where physical properties do not depend 
on the initial preparation of the system. In the strong 
dephasing limit, Eq. (pG) is reduced to the simple form 



TrGoo((9,r) ::^exp 



:rrr(l - z') 



(22) 




:t]:n5 



FIG. 6: Comparison among transient-state counting statistics 
Pt=o(N,r) at FrT = 10 in the strong dephasing regime with 
Fr/F = 0.1 for the three initial conditions a) aoo = l,crii = 
(722 = (dashed line), b) an = 1, croo = (J 22 = (dotted line) 
c) (722 = l,croo — (Jii =0 (continuous line) . 



It gives the probability distribution function for the 
transmitted charges 



Poo(2A/- + l,r) =0 
Poo(2A^,r) = ^^^^^^^ exp 



FrT 

2 



(23a) 
(23b) 



Poo{J^) shows a strong even-odd asymmetry: for even 
the distribution is Poissonian, but the probability that an 
odd number of electrons has passed is negligible. Below 
we will see that this strong parity effect manifests itself as 
an enhancement of zero-frequency shot noise. We leave 
the physical interpretation of the parity effect until we 
discuss shot noise in Section 0. 

In the weak dephasing limit, Eq. (pTI) is reduced to 



TrGoo(^,r) =exp — 



cosh 
so that 
Poo(2A^,r) =exp 

Poo(27V-l,r) = 
where 



3rr 

4 

rr/(z) 



1 



3/(^) 



^%inhIlM 



3Fr\ 



4 d 



(24) 



P7v(r), (25a) 



exp 



3Fr \ 



27ri 



1^1=1 



dz 



sinh 



(25b) 



(26) 



This distribution function shows a much weaker (but still 
finite) even-odd asymmetry than the previous case [cf. 
Eq. (p3|)]. Furthermore, the distribution clearly deviates 
from a Poissonian function, indicating that the presence 
of the strong coherent tunneling of Cooper pairs tends to 
correlate the quasi-particle tunneling events across the 
right junction. This is further reflected in the deviations 
of the current noise from the classical shot noise value 
(see discussions in Section ^). 

One may expect that for a long waiting time (Fr 
00), implying very large numbers of tunneled charges, 
the distribution of N should approach a Gaussian. In 
particular, this becomes an exact result if the distribution 
is Poissonian. In our case, on the other hand, we have 



Poo(27V,r^oo)^-PG(7V,r), 



Poo(2A/--l,r^oo) 



where Pq is a Gaussian distribution. 



Pg{N,t) = 



1 



exp 



{N-Ir/2ey 
2r]r 



(27a) 
(27b) 

(28) 



with 77 = 20/27. The distributions for both even and 
odd N are separately Gaussian, but Poo{N,r F~^) as 
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FIG. 7: Stationary-state distribution Poo{N,r) at Fr = 100 
(a) Ti = = r = IQ-^Ej and (b) Ti = Ts = T = 30Ej. For 
a comparison, a normal distribution function given in Eq. ( |2^ ) 
is also shown (dashed line). 




FIG. 8: Counting probability Poo{N, 
\f2Ej (filled circle). For a comparison 
bution is also plotted (empty circle). 



r) at Fr = 4 for F 
the Poissonian distri- 



a whole is not, since even-odd asymmetry is still present. 
In Fig. R we compare the stationary-state results for 
Poo ^ in the weak and strong dephasing limits with 
the Gaussian distribution. 

Finally, it is interesting to understand what happens to 
the stationary counting probability Poo{N^r) in the in- 
termediate regime, i.e., when the dephasing rate is com- 
parable to the Josephson energy. Unfortunately, an ana- 
lytic expression is not available in this case; the numerical 
results, however, are shown in Fig.||, where one can see 
that the distribution function deviates significantly from 
a Poissonian, being suppressed (enhanced) for odd (even) 
N. 

V. SHOT NOISE 

A deeper insight into the transport process can be ob- 
tained in the frequency domain, from a careful analysis 



of the spectral power of current fluctuations. The zero- 
frequency shot noise could be directly determined by the 
second moment of the counting probability Eq. (18), see 
Ref. ||. Here, however, we follow a different route which 
allows us to get the entire current spectrum. To this end, 
we define the noise spectrum as 



S{u;) = \im I dr e^^^ {{61{t + r), 51{t)}) , (29) 

where SI{t) = I{t) - {I{t)) and {A, B} = AB^BA, The 
total current I{t) through the system is related to the 
tunnelinn currents Il/r = —edtfiL/ji across each junc- 
tion byll2l 



m = ^iR{t)-^iut). 



(30) 



To simplify the evaluation of 5'(a;), it is convenient to 
introduce the spectral densities of currents flowing across 
the individual junctions and the cross correlation spectral 
powers. Therefore, in a way analogous to Eq. (P9| ), we 
write (z, j = L^R) 

/CO 
dre'"^^ {{Shit + T),SIjit)}), (31) 
-00 

which allows to express the total shot noise spectrum in 
the form 



(J2 



C 



\SLR{uj)-rSRL{uj)\. (32) 



In the stationary state (/) = (/l) = — (^i?), so that 
Sibj) = Sll{^) = Srr{u;) in the zero- frequency limit. In 
the opposite liftMr^^ ^ 00), S{uj) = {Cl / C^) S rr{uj) = 
(C£/C|) 2e (/paO. In our case, the left junction is 
(nearly) at resonance for the Cooper pair tunneling and 
therefore limc^^oo = 0. 

In order to obtain Sij (uj) we have to calculate two-time 
correlation functions. We follow the standaxd procedure 
based on the quantum regression theoremE£l and define 
the auxiliary matrices 



qp 



7/^^)(t,r) = Tr {e-^^>totWn,e^''"} 



(33a) 
(33b) 



where the index j runs over the left and right junctions 



(j = L^R). These auxiliary operators, x ^ind 77 



satisfy a master equation of exactly the same form as 
Eq. (0) (but with respect to r instead of t and, of course, 
with different initial conditions). Their relevance can be 
understood by noticing that the correlation functions pan 
be expressed directly in terms of their average valuefj: 
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{{SIji{t + T),dlRit)})=e''{dr-dt)Y^ ^ r„ {n,nR\/''\t,T)+rj<^''\t,T)\n,nR) 

riR n=l,2 

- 2{lR{t)f-2e{lRit))5{T), (34) 
{{5lL{t + t), 5lL{t)}) = 2e\dt - dr){dr + T) ^ (2, n«| x'^-^t, r) + r,(^)(i, r) |2, n«) - 2 (/^(i))' (35) 



and 



({<^/L(t + r), (5/fl(t)}) + {{,5/fl(t + r), (^/^(t)}) = e2((9t - dr) Y.W^r + T) (2, nfl| x^^^ + 



|2,n^)- ^ r(n,n^|x^''^+r/(^M^,^i?)|-4(/LW)(/i?(t)) . (36) 

n=l,2 J 



The problem is now reduced to (a) solving a master equations of the form given in Eq. (|^) to get x^'^H^^r), 
and ri^^\t^T) for respective initial conditions, and (b) evaluating Eqs. (jsj), (|3^), and ( [36| ) to obtain the correlation 
functions. Following this procedure and performing the Fourier transforms of the resulting correlation functions, we 
find that in the stationary state 



2e(/) 



^2 + M2 



(37) 



^"^"^-2(i.|-i^(M|C)-2r|A)). 



2e{I) 



and 



2e{I) 



;2 + M2 



^2 + M2 " 



(38) 



(39) 



where we have used the bra-ket notations 



M 



F F iEj/2 -iEj/2^ 

F -iEj/2 iEj/2 

iEj/2 -iEj/2 T/2^ie 

\-iEj/2 iEj/2 T/2-i€j 







(l 


PI 


,\B) = 




[v 












and 



r 



Ej 



/O 



-i 
\+iJ 



(40) 



(41) 



(42) 



A. Zero- Frequency Noise 



From Eqs. (g7D, (|38|), (|39|), and (|32|), it follows that the 
zero-frequency noise is given by 



S{0) 
2e {!) 



-4£2)2 



(43) 



In the strong dephasing limit (F Ej), the second 
term in Eq. ( p3| ) becomes negligibly small, as it vanishes 
as {Ej/r)'^. Therefore, the zero-frequency shot noise is 
enhanced approximately by a factor 2 compared with its 
classical value, 2e (I). This can be underst4iad4n terms of 
the Josephson quasi-particle (JQP) cycletZlt§ll3. Because 
of the fast quasi-particle tunneling across the right junc- 
tion, each Cooper pair that has tunneled into the central 
island breaks up immediately into quasi-particles, and 
quickly tunnels out. The charge is therefore transferred 
in units of 2e (compared with e in classical charge trans- 
fer) for each JQP cycle. This was already confirmed in 
the counting statistics of the transmitted charges. Ac- 
cording to Eq. (|23), the probability that an odd number 
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of electrons are transferred is zero and charges are trans- 
ferred only in pairs. In the weak (F <^ Ej) and moderate 
(r :^ dephasing limits, the semiclassical JQP picture 
breaks down and we do not have shot noise enhancement 
any longer. 

In the limit T <C Ej^ the period of oscillations of 
Cooper pair is very short compared to the typical time for 
for quasi-particles to tunnel out of the central island. The 
system can be regarded as a single-junction circuit, where 
the quasi-particle tunneling events are independent. The 
Fano factor 6'(0)/2e (/) ~ 10/9, becomes much closer to 
unity in this limit. The small deviation from the Poisson 
value is due to the fact that the quasi-particle tunnel- 
ing events cannot be considered as independent because 
Coulomb blockade allows only one Cooper pair to oscil- 
late coherently across the left junction. Therefore, the 
tunneling process corresponding to |2) |1) is likely to 
be followed by |1) |0). It is clear that this behavior is 
related to the residual even-odd asymmetry we found in 
the counting statistics. Eg . (^ ), even in the long waiting- 
time limit (Fr > 1), Eq. (|27^ 

With moderate dephasing (F :^ Ej)^ quasi-particle 
tunneling events across the right junction are strongly af- 
fected by the coherent oscillation of Cooper pairs across 
the left junction. Indeed, this effect gives rise to the 
significant deviation from the Poissonian distribution of 
the tunneling statistics, Eq. (25). Most remarkably, it 
leads to a suppression of the shot noise. The strongest 
suppression, by a factor of 2/5, is achieved at resonance 
{e = 0) for F = \/2Ej^ see Fig. ^ This is reminiscent 
of the shot noise suppression in (non-superconducting) 
double-junction systemsD, whose maximal suppression is 
by factor 1/2 for the symmetric junctions. We empha- 
size, however, that in the latter case, the coherence was 
not essential. In our case, on the contrary, the role of co- 
herence becomes evident by noticing that the dip in Fano 
factor disappears when moving away from the resonant 
condition as shown in Fig. 0. 



B. Finite-Frequency Noise 

In Figs. |l^ and [l^ we show the typical behavior of the 
finite-frequency noise spectrum in the (a) strong and (b) 
weak dephasing limits. It is interesting to notice that 
(only) in the weak dephasing limit (F <C Ej)^ there is a 
resonance peak of the form 



2el 



2C| {uj ■ 



2e' 



■c.o)2 + r74' 
where the resonance frequency is given by 



r2(E2 ^3^2/2) 



4(£:2 +£2)2 



(44) 



(45) 



Clearly, the peak is an effect of coherent quantum os- 
cillations between the two energy levels separated by 

Uq ~^ > Y^£^j + , induced by the Josephson effect 




FIG. 9: Normalized zero-frequency shot noise for e/Ej — 
0, 0.25, . . . ,5. The dip in the noise is most pronounced at 
resonance (e = 0). 





FIG. 10: Current noise power, (a) at resonance {e = 0) for 
r/£^j = 0, .25, • • • , 5, and (b) at a fixed weak dephasing rate 
(r = O.l^j) for e/Ej = 0, .25, • • • , 5. 



across the left junction. As expected, the resonance peak 
is reduced in its height and broadened in its width with 
increasing F. On the contrary, as e increases, the peak 
gets sharper and the peak height increases quadratically 
with £. However, this should not be confused with the 
zero-frequency case, where £ effectively enhances the de- 
coherence effects. As e increases, the Josephson oscilla- 
tion across the left junction becomes faster, and there are 




2 3 



FIG. 11: Typical behavior of noise power spectrum 5'(u;) as a 
function of frequency uo in the (a) strong (ri,2 ^ Ej) and (b) 
weak (ri,2 ^ Ej) quasi-particle tunneling limits. For both 
plots, £ = and Cl — Cr — were assumed. 
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less chances that it is interrupted by the quasi-particle 
tunnehng across the right junction. This, in turns, im- 
phes that the coherent oscihation is better defined and 
the spectral component (especially Sll{^)) at frequency 
ooq is highly enhanced. For a vanishingly small quasi- 
particle tunneling rate, ^ll (^)/2e/ would approximately 
become a delta-like function, centered at = ^/EJ -\- e"^. 
However, one should not be misled by this result, since 
the noise is always proportional to the average current, 
which vanishes in this limit. 

It is worth mentioning here on the relation between 
this result andj-die description of the noise output from 
linear detectort2l. In the setup considered in this work, 
the right electrode has the role of the detector of Cooper 
pair oscillations; since the total current in the circuit is 
due to quasi-particle tunneling (i.e. it is a dissipative cur- 
rent), the output signal may be considered as classical. 
This has to be compared to the case of a detector mea- 
suring the charge on the island. There, the back-action 
of the detector was essential to produce an observable 
result. In the case of the current, instead, the "detector" 
is intrinsically part of the system and it couples to the 
observed quantity in an essentially non-linear way. 

VI. CONCLUSIONS 

In this paper we considered properties of the distribu- 
tion of the transmitted charge in a superconducting SET 
tuned close to a Cooper pair resonance. The dominant 
process to the transport in the regime considered here, 
is the JQP cycle, a process in which coherent Cooper 
pair oscillations are accompanied by (incoherent) quasi- 
particle tunneling. The interplay between the coherence 
and the strong Coulomb blockade manifests itself in var- 
ious ways both in the counting statistics and in the shot 
noise. We found two distinct regimes characterized by 



different ratios of the time scales for dephasing and re- 
laxation, Tr <C r^^ in the strong dephasing limit or ~ r^p 
in the opposite case of weak dephasing. 

A generic feature of the counting statistics, valid in 
both the regimes, is its even-odd asymmetry related to 
the fact that charge transport is mediated by the Cooper 
pair tunneling. Other properties are more pronounced in 
one of the two regimes. An example is the dependence of 
Pt{N^ t) on the initial time t. This is clearly visible in the 
strong dephasing limit while quickly lost in the weak de- 
phasing regime since, because of the strong Josephson en- 
ergy, the state changes significantly before quasi-particles 
have any time to be produced. Another important point 
is that the counting statistics is not Poissonian, due to 
the relevance of correlations between different tunneling 
events. As a consequence the Fano factor is different from 
the classical value. The maximal suppression of the zero- 
frequency shot noise is observed when the quasi-particle 
tunneling rate is comparable to the frequency scale of the 
coherent Cooper pair oscillations. 

We finally investigated the shot noise at finite frequen- 
cies, which shows a resonance peak at the Josephson os- 
cillation frequency. This maximum can be interpreted 
as an effect of coherent quantum transitions between the 
two energy levels involved in the transport phenomena in 
the device. 
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